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ABSTRACT 


The effects of solar radiation pressure as the main environmental 
disturbance torque have been incorporated into the model of the rigid 
orbiting shallow shell and computer simulation results Indicate that 
within the linear range the rigid modal amplitudes are excited in 
proportion to the area to mass ratio. The effect of higher order terms 
in the gravity-gradient torque expressions previously neglected has 
been evaluated and found to be negligible for the size structures 
under consideration. A graph theory approach has been employed for 
calculating the eigenvalues of a large flexible system by reducing 
the system (stiffness) matrix to lower ordered submatrices. The 
related reachability matrix and term rank concepts are used to verify 
controllability and can be more effective than the alternate numerical 
rank tests. Control laws have been developed for the shape and orienta- 
tion control of the orbiting flexible shallow shell and numerical results 
presented. Hybrid control systems consisting of point actuators and 
passive articulated dumbbell damping devices are evaluated for possible 
use with systems nominally oriented along the local horizontal and can 
result in improved fuel consumption and transient characteristics as 
compared with the active control system operating alone. 
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I. INTRODUCTION 


The present grant represents a continuation of the effort attempted 
in the previous grant years (May 1977-May 1980) and reported in Refs. 

1.6.* Attitude control techniques for the pointing and stabilization 
of very large. Inherently flexible spacecraft systems are being inves- 
tigated. Control of both the orientation relative to the local vertical - 
local horizontal frame of reference and also the shape of the surface 
will often be required to meet the objectives of the various proposed 
LSST missions. These include: communications systems; earth observation 

and resource detection systems; orbitally based electronic mail systems; 
and also orbitally based platforms to be used for the collection of 
solar energy. 

Figure 1.1 illustrates a plan of development of a system software 
capability for use in the analysis of LSST dynamics analysis. This plan 
can be subdivided into four different components: (1) system dyanmics; 

(2) structural dynamics; (3) application of control algorithms; and 
(4) the simulation of enviornmental disturbances. 

Within the system dynamics block, a general formulation of the equation 
of motion of an arbitrary flexible body in orbit has been developed and 
includes also the effects of the Earth's gravity-gradient within the model. 
From this the equations of motion have been developed for homogeneous beams, 
plates, and shallow spherical shell-type structures. 

*For references cited in this report, please see list of references 
after each chapter. 


1.1 






Fig. 1.1: DEVELOPIIEOT OF SYSTZH SOzTvJASZ FOR 

LSST DYI^AMICS ANALYSIS 



/ operational 
IP- in progress 
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For the simpler homogeneous structures such as unconstrained beams and 
circular plates the mode shape functions and elastic frequencies are 
available as relatively simple analytic expressions. For more complex 
systems such as a square or rectangular plate or non-isotropic systems 
these expressions must be calculated based on numerical techniques, usually 
finite element algorithms or numerical computations of rather complex 
analytical approximate relationships. For the complete simulation of the 
controlled dynamics of such systems, the frequencies and the modal shape 
functions of the retained modes are necessary inputs. 

Once the uncontrolled dynamics and stability of these systems has 
been investigated various control algrithms are then used to establish 
controllability and develop control laws for the number and location of 
actuators assumed in the model. During the last grant year the acquisition 
of the ORACLS^ software package, developed by NASA-Langley , has proved to 
be extremely useful in enabling us to develop control laws based on pole 
placement, decoupling, and optimal control theory techniques. 

In accordance with our proposal to NASA dated January 15, 1980®, we 
have also initiated a preliminary study of the disturbance forces which are 
primarily attributed to the effects of solar pressure forces and moments 
and the related problem of thermally induced deflections. The items 
indicated by a check mark have been essentially completed by the end of 
this grant year while those indicated by (IP) are still in progress at 
this time. 
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Chapter II of this report describes the results obtained to date 
in the study of the environmental disturbances and also an examination 
of the effect of higher order terms in the development of the gravity- 
gradient torque expressions which were previously neglected, in response 
to a question asked during last year's final oral presentation. 

A paper presented at the Third VPI&SU/AIAA Symposium on the Dynamics 
and Control of Large Flexible Spacecraft, June 1981, forms the basis 
of Chapter III and illustrates how graph theoretic techniques can be 
used to facilitate the calculation of eigenvalues for large order systems 
and also to reduce the number of computations required to verify system 
controllability . 

In Chapter IV control laws are developed for the shape and orientation 
control of the orbiting flexible shallow shell and typical required con- 
trol forces are simulated. Attention is also given to the calculation 
of the modal frequencies of the shell which are seen to depend on the shell 
curvature as well as the material properties. 

A second paper to be presented at the 1981 AAS/AIAA Astrodynamics 
Conference, August 1981, forms the basis of Chapter V, and examines the 
use of combinations of point actuators together with passive articulated 
dumbbell damping devices for large systems nominally pointing towards 
the Earth (i.e. - oriented along the local horizonal) . The objective 
here is to demonstrate the improvement in fuel consumption and transient 
response characteristics which could be realized with such hybrid control 
systems. 

Chapter VI describes the main general conclusions together with 
recommendations for further work. 
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Chapter II 

ENVIRONMENTAL DISTURBANCES 

A spacecraft in orbit is subjected to many disturbance forces and 
torques which arise from the interaction between the spacecraft and 
the surrounding environment. The most common sources of such environ- 
mental disturbances are: 

1. the earth's atmosphere; 

2. the earth's magnetic field; 

3. the earth's gravitational field; 

4. meteorite impacts; 

5. earth shine radiation; 

6. solar radiation; and 

7. influence of other planets and the sun. 

At the orbits in which the proposed future space missions involving 
large space structures are flown, it is expected that the major source of 
disturbance is the solar radiation. However, since the dimensions of the 
proposed large structures are an order of magnitude larger than those of 
more conventional spacecraft, a reinvestigation into the gravity-gradient 
force and torque expressions which include higher order terms and higher 
harmonics of the earth's gravitational field is needed. Further, if gas 
jets are used for attitude and shape control, unwanted forces and torques 
may arise from the internal leaks in the system. 

In section 2.1, the effects of point forces on a long, flexible beam 
are discussed. The method of approach can be extended to other types 
of structures. Also, some of the remarks made here also hold true, in 
general, for other types of structures. 
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An order of magnitude analysis of the gravity-gradient forces and 
torques due to the previously neglected higher order terms and higher 
harmonics of the earth's gravitational field is presented in section 
2.2 for the class of large space structures currently proposed. 

The effects of solar radiation pressure on the types of structures 
considered in Refs. 2. 1-2. 3 (i.e. beams, plates and shallow spherical 

shells) are discussed in section 2.3. 


2.1. Point Force Acting on a Long, Flexible Beam 

Consider the case of a long flexible orbiting beam which is nominally 
oriented along the local vertical, as shown in Fig. 2.1. The linearized 
equations of motion of this system with disturbances acting were obtained 
in Ref. 2.1 as. 


0 " + 30 = C / J 0)^ 
y y c 


( 2 . 1 ) 


e" + (“ /“ = E /M (n = 1,2,...) (2.2) 

n ncn nnc >»/ 

where, 9 = pitch 

e = nondimens ional modal amplitude (A h) 

n n 

= beam natural frequency; co^ = orbital angular velocity 
= length of the beam; = n^^ modal mass of the beam 
Cy = pitch disturbance torque 

^n ~ modal component of the external disturbance forces 

Jy = pitch moment of inertia of the beam 
( )'= d/dx and T=u^t 

t = time 

It is assumed that a point force of magnitude, p 2 ^(t) is acting on the 
beam at ? = as shown in Fig. 2.2. 
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The torque and the modal force components in eqns. (2.1) and (2.2) are given by, 


where, (^^^^( 5 ) = n*-^ mode shape function of the beam. 

By recalling, J^, = m£^/12; = m, where m is the mass of the beam. 


one arrives at. 


Wi ' ' Wc*- 


where 


, fjL = F^/mu)2jl ; ? =C^/Jl + H 


Thus, one can rewrite eqns. (2.1) and (2.2) as, 
0 " + 30 = 

e" + (o) /(D )^e = <{)^°^f, 

n ' n c n 1 


(2.3) 


(2.4) 


where, 

The system of eqns. (2.3) and (2.4) may be represented by equivalent 
spring mass systems with forcing functions. Hence, one can define the 
following quantities for the above system assuming f^ = a^^sin ( ^ t) , 
where, m = frequency of the point force, and x = m^t 

6 g = 4(?-!5)a^ ; = (}i^^^a^/(a)^/io^)^ (2.5) 

Q 0 g and represent the deflections in each of the component modes 
of the system when a constant force of f^^= a^ is impressed on the systemTj 

By using eqn. (2.5) one can find the relative magnitudes of the influence 


of on each of the elastic modes for a given position of F^. Fig. 2.2 

shows the values of e plotted against the position of F, for the first four 

iis 1 

elastic modes of the beam. 
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From this figure one can conclude that if the force is located at the 
nodal point of a mode, then that particular mode is not influenced by F^. 
Also, one notes that, except for small regions near the nodal points of 
the first and second modes , for any postion of F^ along the beam, the first 
and second modes are influenced quite significantly as compared to the other 
modes. Similar diagrams can be drawn when more than one force is acting 
on the beam. In the design of a controller which uses point actuators 
on the structures, such diagrams may be useful in a preliminary determination 
of actuator locations. 

2.2. Gravity-Gradient Forces 

In the studies presented in Refs. 2.1-2. 3, the force at any point in the 

body due to the gravity-gradient was expressed as, 

1 = 1 + M r (2.6) 

o 

where, f = gravity force intensity at the given point (force/mass) expressed 
in the body principal axes frame; 

f^ = gravity force intensity at the center of mass of the body expressed 
in the body frame (force/mass) ; 

r = instantaneous position vector of the given point with respect to 
the center of mass and expressed in the body frame; 

M = matrix operator defined in eqn. (A. 10) (Appendix A) 

A brief development of eqn. (2.6) is presented in Appendix A. One can 
note that the terms containing (Ap/p)^, (Au)^ and (Ari)^ and higher powers 
are neglected in the derivation of eqn. (2.6) based on the assumption that 
they are very small as compared to unity. 
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One can estimate the order of magnitude of these quantities by noting that, 

(Ap/p)ro ) ; An 0'(r^); Am-'^ , where, I is the largest charac- 

o 

teristic dimension of the structure and, p^ is the orbital radius. For 
currently proposed future large space structures the largest characteristic 
dimension is in the range of 100m to 10km. Hence, even if we consider the 
structure to be in a low altitude earth orbit (e.g. pj^7000km) , the magnitudes 
of (Ap/p), An and Au will be in the order of 10~^ to 10~^. Thus, the 
quantities (Ap/p), An and Am are still much less than unity. Hence, one can 
neglect terms with (Ap/p)^, (Am)^, (An)^ and higher powers in the gravity 
force expression without introducing any significant error in the mathematical 
model . 

The matrix M in eqn. (2.6) is the matrix defined in eqn. (A.9) and is 
given by. 


M= M^°) + ? ^ < 1 

/o) = 


(2.7) 


where, ■' = matrix defined in eqn. (A. 11) of Appendix A; and 


M 


(s) _ v_a^ „ „ t, X 

- p3 B (T^T2) 


-1 


o 

The term corresponds to a spherically symmetric gravitational field. 

The contributions due to higher harmonics of the gravitational field are 
Cs) 

represented by, M (s = 1,2,...). 

If we consider the contribution due to the zonal harmonics, one can 
easily show that, 0 if the origin of the inertial coordinate frame 

is chosen at the center of mass of the earth. The contribution due to the 
second zonal harmonic is given by. 



(T^T2) 


-1 


( 2 . 8 ) 
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where, 


3(3c2n+l) 

6c2n 


0 

0 


B 



6c2ri 

- 1 (7c2n+l) 

0 0 -3(c2n+ i ) 

2 


and J 2 = 1.08263 x 10 , the second zonal harmonic coefficient. 

A comparison of the elements of with the elements of show 

( 2 ) 

that each element of M is three orders of magnitude smaller than the 
corresponding element in Further, it has been estimated by Bainum^'^ 

that the percent error involved in evaluating force or moment expressions 
by neglecting the earth's oblateness (i.e. second zonal harmonic) is about 

_3 

3.76 X 10 at synchronous altitudes. Thus, one can conclude that the 
contributions of higher order terms and higher harmonic components of 
the earth's gravity potential in the gravity-gradient force and torque 
expressions are negligible, even for the currently proposed large space 
structures. 


2.3 Solar Radiation 

At the orbital altitudes suggested for the proposed future large 

space systems one of the major sources of environmental disturbances 

2 5 

would be due to solar radiation pressure. Karymov * has developed 
vector expressions for the forces and moments due to solar radiation 
pressure on a body. It is shown that the magnitude of the force and of 
the moment due to solar radiation on a sphere is independent of the direction 
of the incident radiation and also the surface properties of the sphere. 
Several investigators have considered the effect of solar radiation pressure 
on the dynamics of spacecraft?’^ 2.12 


2.6 



The majority of the spacecraft considered by these investigators consisted 
of a small central rigid satellite to which flat plate appendages were 
assumed to be attached. Since, the centers of the solar radiation pressure 
on these appendages are not usually located at the center of mass of the 
entire system, torques are exerted on the system. One can utilize the 
moments thus generated for attitude control of the spacecraft,^ * ^ 2.10 
controlling the orientation of the plates and/or vanes which can rotate at 
the ends of the appendages. In the case of large space structures, due 
to the Inherently large surface area of the system, a significant amount 
of disturbance torque may be generated due to solar radiation pressure. 

In addition, the deformations induced by thermal stress due to solar heating 
may give rise to additional disturbance torques due to the interaction of 
the incident solar radiation on the deformed surface. A preliminary analysis 
of these effects of solar radiation on a rigid orbiting shallow spherical 
shell is presented here. To the authors' knowledge this is the first such 
analysis to consider the effect of solar radiation on an orbiting shallow 
spherical shell. 

The expression for forces and moments due to solar radiation pressure 
on a shallow spherical shell structure are derived using the approach pre- 
sented in Ref. 2.5. One can easily show that the moment due to solar radiation 

pressure on rigid beams and flat plate structures with uniform surface pro- 

2.5 

parties is zero. The force due to solar radiation pressure is given by, 

? = / dt^ (2.9) 

s 

— oh 

df^ = fgtifg ; df^ = w(l+e'^+2ecos2v) cos v ds 


where , 



e = reflectivity 

T = unit vector in the direction opposite to that of incident light flux 

A 

n = unit outward normal to the surface 

7T A A 

•j “ V = angle between x and n 

w = light energy density = 4.72x10 gms/cm^ (for spacecraft orbiting 

in the vicinity of the earth, i.e. at or near the earth’s dis- 

tance from the sun) . 

The forces and torques due to solar radiation pressure on a shallow 
spherical shell are derived in the following section. 


2.3A. Moments Due to Solar Radiation Pressure on a Rigid Shallow 
Spherical Shell 

Fig. 2.3 shows a shallow spherical shell with the direction of incident 

radiation indicated by -x. The moment due to solar radiation pressure can 
2 5 

be expressed as, ' 


M = (l-e)M^ + eM^ 


( 2 . 10 ) 


s ' ' s s 

where, = moment as calculated by assuming a totalling absorbing surface 
i.e. e = 0; = moment as calculated by assuming a totally reflective 

surface, i.e. e=l. 


Evaluation of : 
s 


given by. 


The moment, for a completely illuminated shallow shell (Fig. 2. 3) is 
2.5 


= ho S2 r^(Txi') (2.11) 

—8 

where, h^ = constant (=4.64x10 N/sq. m) 

S 2 = projected area of the shell on a plane normal to x 

i" = unit vector perpendicular to x in the plane of i and x 

Xq = distance between the center of mass and the center of pressure. 
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that. 


After expressing, x = a^i+b^j+c^k, one can easily derive from Fig. 2.3 


and. 


S2 = = I 

(xxi ) = (l-a 2 )'^(c^j-b^k) 


( 2 . 12 ) 


After substituting eqn. (2.12) into eqn. (2.11) one obtains, 


= iTi^h ^ a (c j-b k) 
s o 2 o' o'^ o 


(2.13) 


Evaluation of : 
s 


The moment due to solar radiation pressure on a totally reflecting 
surface is given by,^'^ 

= 2h /(nxr) (x.n) ^ds (2.14) 

s o 

« s 

where, n = outward normal to the surface of the shell at any point: 
r = position vector from the center of mass of the shell. 

One can express the vectors in eqn. (2.14) as (see Fig. 2.4), 

A /V ^ A 

r = - pi+Re ; n = e (2.15) 

n n 

A A A 

X = (a c6,+b„c6sp,+c^sBs3,)e +(-a„sE. 
olo 10 In 01 

A A 

+b^ce^c3+c^s3c3^)eg +(-b^s3+c^c3)6g 

A AAA 

where, p = R-H/R ; e^ = c3j^i+c3s3j^J+s3s32^ 

A A A A A A A 

eg =-Sg iic3c3j_j+s3c3^k ; e^ =-s3j+c3k 
1 .. 

a ,b ,c = direction cosines of x in the body frame, x, 

000 , 

Hence, " — ^ o '' 

nxr = Ps 3 j^eg 

A A 

(x*n)2 = (a^c3j^+b^c3s3jL+CQs3s3j^)^ J (2.17) 

After substituting the results of eqn. (2.17) into eqn. (2.14) one 
obtains, „ _ „ 


(2.16) 


6 ' 

1 

27T 

(a c 

/ 

/ 

33^=0 

3=0 

\J 

+c 

s3s3,)^ e 


where, = maximum value of 3j^ (see Fig. 2.4) 


2.9 



(2.18) 


After performing the required integration, 

2 

? - uh„R3(l-H/2« [ 2 I - |h 

Eqn. (2.18) can be further simplified for small values of HA by 
noting that, H/R = (as R =Jl^/2H see Fig. 2.5), i.e. 

= IbTTh r 3 ( Y )^a (b k-c j) (2.19) 

s o il o oo 

After substituting the results of eqns. (2.13) and (2.19) into eqn. 
(2.10) and noting, R/ £ = £/2H, (as R =£^/2H, see Fig. 2.5), one obtains 

Mg = Y h^£3(H/£)(l-5E)a^(c^j-b^k) (H/£«l) (2.20) 

If the concave surface of the shell is facing towards the Sun, then 
M^ is obtained by multiplying the expression given in eqn. (2.19) by -1. 
Thus, for the case of the sun shining on the concave surface of the shell 

Mg = Y h^£3(H/£)(l-3e)a^(c^j-b^k) (H/£«l) (2.21) 

If the shell is moving in the earth's shadow, then, M^=0. Fig. 2.6 
shows the geometry of the earth's shadow. For low altitude earth orbits, 
sufficiently accurate results can be obtained by approximating the earth's 
shadow to be a cylindrical shadow of radius equal to that of the earth 
(see Fig. 2.6). The following conditions determine whether the spacecraft 
is in the earth's shadow or not. 

(^•p^) < 0 ; and, [I-(t-Po)^]^ 1 ^ (2.22) 

A 

where, p^ - unit vector along the line joining the earth and spacecraft 
and pointing towards the spacecraft. 

Pg - earth radius. 
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If the spacecraft is moving in a circular orbxt and if the spacecraft 


orbital angular position is measured from the spacecraft positon when the 
spacecraft is in the plane containing x and the orbit normal (and T*p^> 0 
then the above conditions can be expressed as, 

I (1-b^)^ I cosT< 0 ; and (sin^x+b^cos^x)^ — *^e^^o (2.23) 

where, x=m^t and b^ = component of x along the orbit normal. 

The forced response of a rigid shallow shell with gravity stabilizing 
dumbbell due to solar radiation disturbance torques can be obtained by 
solving the following set of equations. The homogeneous form of these 


equations have been developed previously in Ref. 2.3 


0" -3fi 0-c y'-k Y = 

y y y y 

y" (1+c^)y'+ {3+k (1 +c^)} Y+3(l+fi )0 = Q 


(2.24a) 

(2.24b) 


4i" >J;-(1-K2 )<{.' = 0 

(p" +4n Hd-n )ip' = 


(2.24c) 

(2.24d) 


6" +c (l+c-)6'+ {4+k (l+c„)}6+4(l-n )(j) 
z z z z z 

-a-n^)r=Q, 


(2.24e) 


where. 


C^®^= (M^*j)/J 
y s y c 


6 a c (l-5e) , a >0 
s o o ’ o 

6 a c (l-3e) , a <0 
s o o ’ o 

0, when the spacecraft is 

in the Earth shadow. 


C^®^= (M *k)/J 0)2 
z s ' z c 


5^ a b (l-5e) , a >0 
s o o ’ o 

5 a b (l-3e) , a <0 
s o o o 

0, when the spacecraft is 

is in the Earth shadow. 


5^ = 2(r£2/m)(H/)0 (~^) 

c 


Q ,Q^ = torques due to solar radiation on the dumbbell. 
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Figs- 2.7a and b show the nature of the disturbance torques, due 
(s) (s) 

to solar radiation, and , for a spacecraft wxth an area to mass 

ratio of 0.02 m^/kg and orbiting in a 250 n. mile altitude circular orbit, 

with the orbit normal oriented at an angle of 45 degrees to the earth-sun 

line. It can be noted from these figures that these torques are periodic. 

(s) 

The average value of the pitch torque, C^, , over one orbit is noted to be 

(s) 

zero. However, it is seen that the average value of the roll torque, , 
over one orbit is not zero. 

If the dumbbell is assumed to be rigid then the net torques on the dumb- 
bell due to solar radiation pressureare zero, l.e., =0 and =0. 

Figs. 2. 8-2. 9 show the response of the rigid shallow shell-rigid dximb- 
bell system subjected to the torques due to solar radiation pressure. The 
following numerical values for the constants have been assumed in obtaining 
these response plots. 


h^ = 4.64x10 ^N/m^ ; (H/ £) =0.02 ; 2. = 50m ; e = 0.5 
In all these cases Initial conditions are assumed to be zero. In Figs. 2.8 
and 2.9 one can observe the beat phenomenon in the roll and yaw responses 
of the system. This indicates that some of the frequencies of the forcing 
functions are close to at least one of the system natural frequencies. One can 
also note that with an increase in the value of area to mass ratio, the ampli- 
tudes of pitch, roll, yaw, y and 6 responses increase in direct proportion. 
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Effect of the Thermal Deflection of the Stabillzina Boom 

Due to the heating effect of the solar radiation, thermal stresses are 
induced in the structure. Thermal deformations are introduced in the structure 
because of these thermal stresses. It is assumed that the surface of the 
shallow spherical shell is controlled to be within the required accuracy 
limits although the control dynamics is not considered here. It is expected 
that the torque due to solar radiation pressure on the deformed shell is 
almost the same as that corresponding to the torque on an undeformed shell. 
However, since it is expected that there will be no control present to control 
the shape of the stabilizing boom, the boom will undergo thermal deformations 
due to solar heating. Additional torques may arise due to solar radiation 
pressure acting on the thermally deformed boom. In the present work, the 
boom is assumed to be a hollow tube of uniform cross section. The following 
additional simplifying assumptions are made in this preliminary investigation! 

1. The boom cross section is of uniform thickness. 

2. Boom deflections occur only in the plane containing the sunline and 
the boom’s undeflected axis. 

3. The temperature of the boom is independent of the boom's axial 
coordinate and the axial boom deflection is of higher order. 

A. The energy exchange between the deflected boom and the satellite 
surface is negligible. 
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5. The earth radiation input may be neglected at the higher altitudes 
under consideration. 

6. The inherent time lags in the heat transfer process are neglected. 

7. Changes of moments of inertia and center of mass shift due to 
thermal deflections are neglected. 

8. Boom deflections occur such that each differential element on the 
boom has the same constant radius of curvature. 

In Ref. 2. A an expression for the torque due to solar radiation pres- 
sure on a thermally bent pair of booms hinged at the center of the unde- 
formed boom set has been derived based on the assumptions (l)-(8). The 

2.4 

expression for the torque is given as * , 


AN = ^ G'D sin 2B. sin B^ 
s 3 1 i 


(2.25) 


where, 


G' = r,(l+-| ) ; D = eq /2kt 

oad j ao '^ '■ 


29 ,^ = length of the boom 


r, = radius of the tube; t, = thickness of the tube; 
a a 

e = reflectivity ; e = emissivity 
= thermal absorptivity 

B^ = Sun incidence angle between the sun direction and 
the undeflected boom axis (see Fig. 2.10) 


The net torque attributed to solar pressure forces on the boom' s 
2.4 


tip masses is given by 
h S£ 


AN = 
m 


D sin 2 B. 

X 


(2.26) 


where, S is the cross sectional spherical area. 
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If it is further assumed for this example that the Sun is shining in 
the orbit plane, then the net torque on the boom attributed to the solar 
radiation pressure on the thermally deflected boom is given by, = 0 and, 

Q = AN +AN (2.27) 

y s m 

By substituting eqn. (2.27) into eqn. (2.24b), one can solve eqns. 
(2.24a) and (2.24b) to obtain the pitch (6) and dumbbell angle (y) responses 
of the system for this case. In the present analysis a numerical inte- 
gration procedure based on fourth order Runge-Kutta method was used in 
solving eqns. (2.24a) and (2.24b). For the example considered here, the 
following data has been assumed, boom material: beryllium copper. 

= 0.15 ; e = 1.04xl0"^/°F ; k = 36.3779k. cals'm-°F-hr; 
e = 0.85 ; q^= 1464.8 k.cals/hr-m^ 

In addition, r,/t, = 1000 ; r, = 5cm ; tip mass = 0.5 (mass of the shell), 

G a u 

have been assimed. 

From the results of the numerical integration, it is noted that, for 
the set of data assumed in this particular example, the effect of the 
torque due to solar radiation pressure on a thermally deformed boom is 
negligible; thus, the pitch and dumbbell angle (y) response in this case ‘ 
is essentially the same as that shown in Figs. 2.8-2. 9 . 

2.3B. Forces Due to Solar Radiation Pressure on a Rigid Shallow 
Spherical Shell 

The force due to solar radiation pressure on a shallow spherical shell 
can be expressed in a similar fashion to eqn. (2.10), i.e., 

Fg = (1-e) + eF^^ (2.28) 
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where, 


= force 


as calculated by assuming a totally absorbing surface, 


i.e. e=0 j 


= force as calculated by assuming a totally reflective surface, 
i.e. e=l. 


Evaluation of F^ : 
s 

The force, F^ for a completely illuminated shell is given by 
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F^ =-h s„T 
s o 2 


(2.29) 


where, S 2 ~ projected area of the shell on a plane normal to x = irjl^a^ . 


Hence 


, F^ = -rh £^a (a i+b j+c k) 

S 0 0^0 o'* o' 


(2.30) 


“■R 

Evaluation of F : 
s 


"“R 2 5 

The force, F^ for a completely Illuminated shell is given by ' , 

x:R 


F = -2h /n(x*n)^ ds 

s Os'' 


(2.31) 

where, n = outward normal to the surface of the shell at any point. 


By using the results in eqns. (2.15) and (2.16), one has, 
-R O 

F^ = -2h r2 r f (a c3,+b cBsB, 
s o „ „ „ „ o 1 o 1 


B ^=0 B =0 


+c^s6s3^)2e^s6j^dBd3^ (2.32) 


where, 6,' = maximum value of 3 . 


After performing the required integration, one obtains. 


'a ■ - I f)*) + ( I )“) i 

+ 2a^(b^J+c^k) ( I )*] (2.33) 

Eqn. (2.33) can be further simplified for small values of (H/£) by 
noting that, (H/R) = 2(H/Jl)^ and (X./R) = 2(H/Jl), i.e., 

Fg -8Trh^R2(H/Jl)2[a2i+2a^(b^j+c^k) (H/t)^] (H/Jl«l) (2.34) 
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After substituting eqns. (2.30) and (2.34) into eqn. (2.28), 
one arrives at, 

A. ^ 

F % -Trh £^[(l+£)a^i+(l-e)a (b j+c k)] (h/JI<<1) (2.35) 

s o o o o o 



If the concave surface of the shell is facing the sun, then, is 
obtained by multiplying the expression in eqn. (2.34) by -1. Hence for 
the case of the sxin shining on the concave surface of the shell, 

-7Th^i^[(l-3e)a^i+(l-e)a^(b^j+c^k)] (H/£<<1) (2.36) 

The expressions in eqns. (2.35) and (2.36) might be useful in estimating 
the spacecraft orbit perturbations due to solar radiation pressure. 


2.17 



References 


2.1. Bainvun, P.M. , Kxjmar, V.K. , and James, P.K. , "Dynamics and Control of 
Large Space Structures - Part B," Final Report, NASA Grant NSG- 
1414, Dept, of Mechanical Eng., Howard University, May 1978. 

2.2 Bainum, P.M. , James, P.K. , Kumar, V.K., and Krishna, R. , "The 
Dynamics and Control of Large Flexible Space Structures - II, 

Part B: Model Development and Computer Simulation," Final Report, 

NASA Grant NSG-1414, Suppl. 1, Dept, of Mechanical Eng., Howard 
University, June 1979. 

2.3. Bainum, P.M. and Kumar, V.K. , "The Dynamics and Stability of Large 
Flexible Space Structures - III, Part B: The Modelling, Dynamics 

and Stability of Large Earth Pointing Structures," Final Report, 

NASA Grant NSG-1414, Suppl. 2, Dept, of Mechanical Eng. , Howard 
University, Sept. 1980. 

2.4 Bainum P.M. , "On the Motion and Stability of a Multiple Connected 
Gravity Gradient Satellite with Passive Damping," Technical Memo- 
randum, TG-872, Johns Hopkins University, Appl. Phy. Lab., Jan. 

1967, (Ph.D. Thesis, Catholic University of America, 1966). 

2.5. Karymov, A. A., "Determination of Forces and Moments Due to Light 
Pressure Acting on a Body in Cosmic Space," PMM, Vol. 26, No. 5, 

1962, pp. 867-876. 

2.6. Crocker, M.Ct, "Attitude Control of a Sun Pointing Spinning Space- 
craft by Means of Solar Radiation Pressure," J. of Spacecraft and 
Rockets, Vol. 7, No. 4, 1970, pp. 357. 

2.7. Donlis, T.J. and Randall, J.C., "A Solar Vane Actuation System for 
Spacecraft Attitude Control," ASME Mech. Conf . , Lafayette, Oct. 1964, 
paper no. 64. 

2.8. Modi, V.J. and Pande, K.C., "Solar Pressure Control of a Dual Spin 
Satellite," J. of Spacecraft and Rockets, Vol. 10, 1973, pp. 355. 

2.9. Modi, V.J. and Tschann, C., "On the Attitude and Librational Control 
of a Satellite Using Solar Radiation Pressure," presented at XXI lAF 
Congress, Konstanz, West Germany, Oct. 4-10, 1970. 

2.10. Pande, K.C., Davis, M.S. and tfodi, V. J. , "Time Optimal Pitch Control 
of Satellites Using Solar Radiation Pressure," J. of Spacecraft and 
Rockets, Vol. 11, No. 8, 1974, pp. 601-603. 

2J.1. Clancy, T.F. and Mitchell, T.P., "Effects of Radiation Forces 

on the Attitude of an Artificial Earth Satellite," AIAA J. , Vol. 2, 

No. 3, 1964, pp. 517-524. 

2.12. Flanagan, R.S., and Modi, V.J., "Attitude Dynamics of a Gravity Gradient 

Satellite under the Influence of Solar Radiation Pressure," Aeronautical J., 
Vol. 74, Oct. 1970, pp. 835-841. 


2.18 




2.20 



0 0.2 0-4 0.6 0.8 1.0 

Normalised distance along the beam, Q =(x/£+ h) 


Fig. 2.2 : Influnce of a Point Force on the Elastic Modes of a Free-free Beam 
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FiR. 2.3: Shallow Spherical Shell with Incident Solar Radiation 
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Fxg. 2.7a; A Typical Variation of the Disturbance Torque, C' ' , due 
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to Solar Radiation Pressure 
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Fig. 2.8a: Forced Pitch ( 0 ) and Dumbbell Angle (y) Responses of the 
Gravitationally Stabilized Shallow Spherical Shell due to 
Solar Radiation Pressure 
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Fig. 2.8c: Forced Yaw Response of the Gravitationally Stabilized Shallow 
Spherical Shell due to Solar Radiation Pressure 




Fig, 2.9a; Forced Pitch (0) and Dumbbell Angle (y) Response of the 

Gravitationally Stabilized Shallow Spherical Shell due to 
Solar Radiation Pressure- Effect of increase in the area 
to mass ratio 



Roll, (}) (rad 









Fig. 2J.0 ; Boom Deflection Geometry 

(only half of total boom assembly shown) 
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Appendix A 


A brief development of the expression for the gravitational force 
per unit mass, f, is presented here. The material presented here is 
extracted from Santini. 


Kinematic Relations 

Fig. A. 1 shows the various coordinate frames used in the development 
of gravity gradient force expressions. O'XYZ is an earth centered inertial 
reference frame (t^) with O'Z along the earth's spin axis and O'X along the 
ascending node. 0ij^i2i2 is the local intrinsic frame centered at the 

body center of mass, 0, with 01^^ along the local vertical and 0i2 perpendi- 
cular to the plane ZO'O (Fig.AJ). OX Y Z is an orbit fixed reference frame 

o o o 

(T2) centered at the center of mass of the body, 0, with OX^ along the local 
vertical and OY^ along the orbit normal opposite to the orbit angular momentum 
vector. Oxyz defines the principal axes (x^) of the body (not shown) . 

Fig. A. 2 shows the sequence of Euler angle rotations adopted in going 
from the orbit frame, body frame, X2(x,y,z). The body frame, 

X2(x,y,z) and the orbit frame, T2(X^,Y^,Z^) are related to each other by the 
following transformation relation, 

T3 ^ T^T2 (A.l) 


where. 


T 


1 


C({ic0 sij)ci(H-c<j)s9si(/ s^si|)-c<jis0c4’ 

-S({)C0 cij)ciJ;-S())s0si|) c<j)Sij;+S(J)S0ciJ^ 

s0 c0si|^ c0ci|; 


The body angular velocity components (m ,0) ) and Euler angular rates 

X y 2 

• • • 

(Q,<p,ip) are related as follows. 


• • 

0) = 0s(})+i(JC(()c0-a) (s4>ci{;+c<j)s0ST|;) 

X c 

• • 

03 = 0c(fi-'l3S!})c6-a3 (c(fiC3j3-S(j)S0si(;) 

Y • • ^ 

03^ = l{3S0+({)+O3^c6S3f3 


(A. 2) 


*Santxni, P., 'Stability of Flexible Spacecrafts', Acta Astronautica, Vol. 3, 
1976, pp. 685-713. 



Fig. A. 3 shows the sequence of rotations adopted in goxng from the 

intrinsic frame, t,„(i,i_i_) to the orbit frame, t_ (X ,Y ,Z ). The 
10 12 3 ’ 2 o’ o’ o"^ 

intrinsic frame, is related to the orbit frame, by the trans- 

formation relation, (see also Fig. A.l) 


"^2 


(A. 3) 


where , 


To = 


10 0 
0 cx sx 
0 - X cx 


Gravity Force 

One can write the gravitational potenial in the most general form as 


where , 


V (p,n, oj) = 


va 


a sS+1, 


+ va E ( ^ ) X (n,m) 
p s=l p s 


(m) 


(A.4) 


fig(n,u) = E (n)cos(mm+(}>gjj^)] 

m=0 

(n) is the m^*^ associated Legendre function of order s. 
and <{)^^ are constants to be obtained experimentally through 
satellite geodesy techniques. 

The gravitational force per unit mass at any point in the body, f, 
with components in the local intrinsic frame is given by. 


where. 


F(p,ri,a)) = Fe + Fe +Fe 
p p n n 0) 0) 


(e , e , e ) are unit vectors in the local intrinsic frame. 
P n ^ 

„ _ 3V _ va^ , “ .a 

^p - 3^ - - X ^ ^ P ^ " 

F = i ^ 

n p 3q P s=l P ® 

m psn 3oj psri s=l p s 
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where, ( )' = 3( )/3n ; () = 3()/3o) 

The gravity force per unit mass at the generic point, P, in the 
body as expressed in the local intrinsic frame at the center of mass of 
the body and for small values of 6p,<Sn and 5m is given by 

3F 


F = Fq + AFq + 


3(p,n,m) 


Ap 

An 

Am 


+ H.O.T 


(A.5) 


where , 



" 0 

- F 

n 

-F sn ” 

m 


Ap 

II 

0 

F 

p 

-F cn 
m 


An 


0 

0 

F sn+F cn 
p n 


_ Am _ 


Fq = gravitational force per unit mass at the center of mass of the body. 
After some mathematical manipulations, it can be shown that. 


.... 3F 

A ^ 

— 0 3(p,n,u) 



B* 



^ psnAm 


(A. 6) 


where , 


B* = 


va‘ 


,( 0 ) 


,(s). 


+ Z ( - ) B""'] 

s=l P 


II 

o 

PQ 

“ 2 
0 

0 

-1 

1 

o o 


O 

0 

-iJ 


,(s) 


(s+1) (s+2)n 


■(s+2)n^ 


-(s+2)|- 

sn 


-(s+2)n^ 

[n'J-Cs+Dn^] 






sn 


[f2' cot n-s^ (s+l+ 
s s 
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To a first order approximation, 

r Ap ^ 

I pspAm j 

where r is the position vector expressed in the body axes frame, 


= (T^T2) 


(A. 7) 


Thus , 


F = Fq + B*(T2^T2)"^r 


(A. 8) 

The gravitational force per unit mass at the generic pointy P , in the 
body can be expressed in the body principal axes frame, by using the 


transformation, (T^T 2 ) in eqn. (A. 7), i.e.. 


where. 


f = fg + M r 


f = T T F 
0 12 0 


M = T^T2B* (Tj^T2) 


(A.9) 


-1 


If we write, M = Z ( — ) M 

Po 


(s) 


(A. 10) 


then. 


jj(0) ^ x^T2B^°\t^T2)"^ ^ 

Po 

^(s) ^ 

o 


For the sequence of Euler angle rotations shown in Fig. A. 1, 


^(0)^ ^ 


3c^({ic^0-l -3s4>c(j)c^0 3c(j)c0s0 

-3s(})ci})c^0 3s^(J)c^0-l -3s(jic0s0 

3c(J)c0s0 -3s<J)c0s0 3s^0-l 


(A. 11) 


It should be noted that the matrices, B^^\ B^^^ and are 


symmetric. 
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Appendix B 


An approximate expression for, x^ (see Fig. 2.5) is derived here. 
Fig. B.l shows the cross sectional view of a shallow spherical shell 
with the various symbols defined. From plane geometry one notes that 
(see Fig. B.l), 

AF . FB = EF . FJ and CG . GD = EG . GJ 

i.e. r^ = A(2R-A) and = H(2R-H) (B.l) 

c 

If it is assumed that, H/A«l (which also implies that, H/R<<1 and 
A/R<<1) one can write, 

A r2/2R and H £^/2R (B.2) 

u C 

Hence, x = H - A % (£^ - r^)/2R (B.3) 

c c 


kxis of 


symmetry 


C 

/ 



surface of a shallow 
spherical shell 

D 






J 


Fig. B.l: Geometry of the cross section of a 

shallow spherical shell. 
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CHAPTER III 


GRAPH THEORY APPROACH TO THE EIGENVALUE 
PROBLEM OF LARGE SPACE STRUCTURES 


Abstract 

A graph theory approach for reducing a matrix to lower ordered sub- 
matrices can be utilized to find eigenvalues of a large flexible system. 
This approach can reduce the computational effort involved in finding 
the eigenvalues. The reachability matrix and term rank concepts are used 
to verify controllability and can be more effective than numerical rank 
tests of the system controllability matrix. A free-free square plate is 
considered as an example to illustrate these techniques. 


I. Introduction 


Large, flexible orbiting systems have been proposed for possible use 
-in communications, electronic orbital based mail systems, and solar 
energy collection. ^ ^ The size and low weight to area ratio of such sys- 
tems warrant the consideration of flexibility as the main contribution to 
the dynamics and control problem as compared to the inherently rigid na- 
ture of earlier spacecraft systems. For such large flexible systems 
both orientation and surface shape control will often be reqiiired. 

An analysis of the effectiveness of control systems for such mis- 
sions is usually divided into: (i) a dynamic analysis of the uncontrolled 
flexible orbiting system; and (ii) control law selection for particular 
types and locations of actuators and sensors under consideration. 

For a dynamic analysis it is necessary to solve the eigenvalue pro- 
blem in order to determine natural (modal) frequencies and mode shapes. 

The eigenvalue problem for such large space structures may pose enormous 
numerical and computational problems for systems characterized by a state 
vector of large dimension. 

The present paper addresses itself to approach such a problem using 
graph theory. The determination of the eigenvalues for many vibratxon 
problems involves the solution to a large dimensional eigenvalue problem. 
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In the literature there are certain methods to solve these problems 
using the special nature of the matrices involved in specific situations 
such as for rotating (gyroscopic) systems.^’** 

In the present paper the methods available to interpret the deter- 
minant of a matrix graphically® are exploited to reduce the higher di- 
mensional matrix into combinations of smaller dimensional sub-matrices , 
wherever possible. The union of the eigenvalues of the smaller matrices 
including their multiplicity results in the eigenvalues of the original 
matrix. ® The reduction procedures involved employ the Boolean equivalent 
of the original matrices to arrive at the smaller dimensional equivalents 
of the original numerical matrix. As the matrix Boolean operations can 
be performed rapidly with comparably little memory requirements with a 
high speed computer, the reduction analysis has definite advantages as 
compared with the direct solution of the eigenvalue problem using the 
original system matrix. The eigenvalues found from the smaller sub- 
matrices will, in general, be more accurate than those calculated from 
the original matrix. 


II. Basic Definitions and Theorems ® 


The following are some of the definitions and theorems useful in 
understanding graph theory as applied to matrices. These concepts will 
be illustrated in detail in Section III where graph theory is used to 
determine the* eigenvalues of a free-free square plate. 

1. Directed graph of a matrix ; D(M) is a directed graph of the matrix, 
M (of order n) with n points a 2 ,..., a^, in which there is a 
directed line from a^ to aj^ if and only if the entry my of the ma- 
trix is not zero, and, in addition, there is a directed line from 
each point to itself. 

2. A directed graph is strongly connected or more briefly, is strong , 
if, for any two points, a and b, there exists a directed path from 
a to b and a directed path from b to a. 

3. A strong component of a directed graph is a maximal strong sub-graph. 

4. The points of any directed graph are partitioned into mutually exclu- 
sive strong components . However, not every line of a directed graph 
necessarily lies in a strong component. 


Graphical Interpretation of a Determinant ; The following assertions are 
usefxil in interpreting the determinant of a matrix graphically. 

Lemma 1 ; The lines of D(M) , corresponding to the entries of M whose pro- 
duct is a non-zero term of |m| (|m| is the determinant of matrix, M) , 
constitute a collection of disjoint directed cycles within the directed 
graph, D(M). 
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LpTnma 2 : Every line of D(M) which occurs in a directed cycle appears in 

a strong component of D(M) . 

Lernma 3 : Each line of D(M) which does not lie in a strong component of 
D(M) is coordinated with an entry of M which does not occur in any non- 
zero terms of [m|. 

Lemma 4; The product of the determinants of the sub-matrices of M 
associated with the strong components of D(M) gives, up to the sign, 
the value of the determinant of M. 

The proofs of Lemmas 1-4 are provided in Ref. 5. 

The following theorem can be used to test whether a matrix can be 
reducible to sub-matrices or not. 

Theorem ^ ; The matrix M=(mj^j) of order n is irreducible if and only if 
n=l or if n>l, and for any i and j such that 1 <^ i£n, l^^n and ij^ j , 
either m^jT^O or there exist i]_,i 2 ,...» is such.~thiat °^,i 2 '-**“ig 

i.e. there is a path between the i and j nodes in a digraph sense. 

Theorem on Eigenvalues ^ ; 


Let E(M) be the collection of all the eigenvalues of M, including 
all multiple eigenvalues counted separately. With the construction of 
sub-matrices, Mj_, used in Lemma 4, then let E(M^) have the corresponding 
■ meaning for the collection of eigenvalues of the sub-matrix, M^. 

Theorem: Using the above notation for E(M) , the eigenvalues of M are 
given by E(M)=E(Mj^)UE(M 2 ) . . .UE(Mj-) , including mtiltiplicity. In words, 
the set of eigenvalues of the matrix is the union (including multi- 
plicity) of the sets of the eigenvalues of the sub-matrices corresponding 
to the strong components of the directed graph of the matrix. (The proof 
of this theorem is given in Ref. 5.) 


III. Algorithm for Matrix Reduction 

An algorithm is given here to find the bub-matrices corresponding 
to the strong components of a directed graph of a matrix M(order n) . 


Step 1 : Replace the diagonal elements (whether zero or non-zero) by 

unity and replace all non-zero non-diagonal elements by imlty. We ob- 
tain an adjacency matrix. A, of the digraph, D(M) . We obsejrve that 
A is a Boolean matrix - i.e. contains only "zeros" and "ones" as elements. 


Step 2 : Find the matrix A by Boolean multiplication where i is either 

determined by A^=A^"^ (i.e. 1+1=1, 1+0=1, 0+0=0 rules) or i=n-l. This 

Mtrix is called the reachability matrix (or path matrix) and is denoted 
by R* 


Step 3 : Let 82 ^ be the collection of all the indices, i, such that there 

appears the number "1" in both the i^h place of the fxrst row of R and 
of the first column of R. 
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These Indices name the points of D(M) in the same strong component as 
a-^. Delete from the matrix, R, all rows and columns whose indices are 
in the set, but preserve the original row and column indices from R 
in the resulting sub-matrix of R. 

Step 4 ; Repeat step 3 on this resulting sub-matrix, obtaining the 
collection of sets of indices 82 , 33 , ... ,S^. Then, clearly S]_,S 2 » • • • ,3^ 
are pairwise disjoint sets of indices whose union is {l, 2 ,...,nj and 
each set, contains all the indices which belong to the points in 
one strong component of D(M). 

Step 5 ; Let be the n^ by n^^ sub-matrix formed from the rows and 
colimns that are in the set S^. Then by the Theorem on the Eigenvalues 
and Lemma 4, the original matrix, M, is completely reduced to the ma- 
trices, M^, and each is irreducible. Also, the eigenvalues of M are 
the union of the sets of eigenvalues of the sub-matrices Mj[ (including 
multiplicity) . 

The numerical effort involved in reducing the original matrix is 
nominal as the matrices involved are Boolean. The multiplications can 
be done row by column wise using an assembly level language. The me- 
mory involved is also nominal as every element can be stored in a bit 
instead of a word. Thus, this procedure is very attractive for large 
dimensional systems. The above graphical procedure is applied to a 
mathematical model of a free free square plate.®’® 


IV. Example of a Free-Free Square Homogeneous Plate 

The transverse displacement of the free-free square plate may be 
represented by® : 

J N 

W(x,y,t) = n^l 

where (j)j(n), ^re free-free beam functions 

r^(t) is a time dependent function 

r = n+(j-l)N with n = 1,2,...N 
j = 1,2,. ..J 
r = 1,2,...JN 

The dynamic equations of the plate are given by®: 

M g+Kg = qf(t) ( 2 ) 

where M = ml, m = mass of the plate 
T T 

8 “ ^^1’ ^2’’*’ ^ ^ ~ [’?]_> *?2’’*'^ 

With p(x,y) as the spatial distribution of the control force, then. 
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a b 

■ q_ = / / P(x,y)<p.(x)4/ (y) dydx 
^00 J a 

and the general element of the stiffness matrix may be expressed 


K = 
rs 

D 

ab 


b^ 

a^ 

J ^2 n 





+ V 


'^nn' 



'•^nn'^ 






s = 

n'+(j'-l)N 






a 




b 

d^i{j 



II 

•r-t 

H 

a f 
0 


■ 

,dx ; I , = b 
’ nn 

/ 

0 

11 

dy^ 


dy 


a 

a / 
0 

dx 

dx 

rO 

II 

1 

b 

/ 

0 

dnb 

dy 

dy 

dy 


a., 0 are the free-free beam frequencies. v is Poisson's ratio, and 
D^is ?he flexural rigidity. The non-zero elements of the K matrix are 
calculated- for m = 276800 kg, a=b=100m. 

= 0^ = 6.0 ; = 0^ = 4.73004 ; “2 ~ ^2 ~ 

= 0^ = 7.85324. 

For J = N = 4 we have 16 modes for the plate. Fig. 1 illustrates the 
location of the non-zero elements in the K matrix. (See Step 1, 
Algorithm) corresponding to the following actual nimerical values 
(Table I) 

TABLE I 


The First Four Free- Free Beam Frequency Parameters : 


0.0,0.0,4.73004,7.85324 
Non-zero elements of K-matrix: 


K(3,3)=34.919 

K(12,4)=-19.9054 

K(14,6)=15.4658 

K(15,7)=-11.6130 

K(16,8)=96.9152 

K(4,10)=17.5471 

K(9,ll)=-5. 31786 

K(12,12)=326.604 

K(6,14)=15.4658 

K(7,15)=-11.6130 

K(8,16)=96.9152 


K(9,3)=8. 03522 

K(10,4)=17.5471 

K(16,6)=17.8584 

K(6,8)=15.4658 

K(3,9)=8. 03522 

K(10,10)=34.9190 

K(ll,ll)=76.8770 

K(7,13)=17.5471 


K(8,14)=56.1775 


K(13,15 

Ka4,16 


) =-19. 9054 
)=96.9152 


K(ll,3)=-5. 31786 
K(6,6)=13.3938 
K(7,7)=34.9190 
K(8,8)=386.913 


K(4,4)=265.337 

K(8,6)=15.4658 

K(13,7)=17.5471 

K(14,8)=56.1775 


K(9,9)=34.9190 K(ll,9)=-5. 31786 

K(12, 10)=-11. 6130 K(3,ll)=-5. 31786 
K(4 ,12) = -19 . 9054 K(10 , 12)=-11. 6130 
K(13,13>= 265.337 K(15 ,13)=-19 . 9054 


K(14,14)=386.913 K(16,14)=96.9152 
K(15,15)=326.604 K(6 ,16)=17 . 8584 
Ka6,16)=1732.85 


3.5 



Col. 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

Row 1 


r 



■ 

■ 

■ 

■ 

■ 

■ 

■ 

n 

■ 

■ 



2 

■ 

■ 

■ 


■ 



■ 

■ 

■ 

■ 

■ 

■ 

■ 



3 

■ 

■ 


■ 

■ 

■ 



X 

■ 

H 

■ 





4 

■ 

■ 


X 


■ 

■ 

■ 

■ 

X 


X 

■ 

■ 



5 



■ 

■ 

■ 


■ 

■ 

■ 

■ 

■ 

■ 

■ 

■ 



6 



■ 

■ 


X 



■ 

■ 

■ 

■ 

■ 

X 



7 

1 

■ 

■ 


■ 

■ 

X 

■ 

■ 

■ 

■ 


X 


X 


8 

■ 

■ 

■ 



X 


X 

■ 

■ 

■ 

■ 

■ 

X 


X 

9 



X 

■ 

■ 

■ 

■ 


X 


X 


■ 

■ 

i 


10 


■ 

■ 

X 

■ 

■ 

■ 



X 


X 


■ 

1 


11 



X 






X 


X 

■ 

■ 

■ 

■ 


12 


■ 

■ 

X 






X 


X 


■ 

■ 


13 







X 


■ 

■ 

■ 

■ 

X 


X 


14 


■ 

■ 



X 


X 

■ 

■ 

■ 



X 

■ 


15 


■ 

■ 

■ 

■ 

■ 

X 

■ 

■ 

■ 

■ 

■ 

X 


X 


16 


■ 

■ 

■ 


-ij 


X 

■ 

■ 

■ 

■ 

■ 

X 




Fig. 1. Non-zero Element Pattern of the K-Matrix. 
(X indicates non-zero element) 











































From Fig. 1 it can be seen that the isolated nodes on the digraph 
correspond to the indices "1", "2", and ”5", since none of these cor- 
respondingly numbered rows and columns contain non-zero elements. Thus, 
the digraph centered at these nodes consists only of a directed line 
from the node to itself. As an example of a strongly connected component 
of a directed graph, we consider the nodes corresponding to the indices, 
"3", "9", and "11". It can be seen that non-zero elements appear in the 
third, ninth, and eleventh places in both the third row and columns, re- 
spectively. Therefore, the directed graph between any two of these 
nodal points is strongly connected. The remainder of the digraph is 
constructed based on similar reasoning. It can be verified that the 
Boolean equivalent^ of is the same as that for the K matrix itself. 

This matrix is thus the path matrix (reachability matrix) and its di- 
graph is shown in Fig. 2. 

The reachability matrix shown in Fig. 3 is subdivided following Step 
3 of the algorithm. For example, there is a "one" contained in the first 
row and the first column of the reachability matrix. No other non-zero 
element appears in either the first row or the first column. Thus the 
subset, Sj_, has, in this case, only the single index "1". By similar 
reasoning the second subset, S 2 , associated with the "one" contained in 
the second row and the second cplumn of the original matrix (the first 
row and the first colimin of the reduced matrix) also has in this case a 
single index, " 2 ". 

We now consider elements appearing in the i*-^ places of the third 
row and column, respectively, of the 14x14 matrix formed by excluding Sj_ 
and S 2 . It can be seen that the nvanber "one" appears in the third, 
ninth, and eleventh places of both the third row and column of the reach- 
ability matrix. Thus the subset, S 3 , now contains the three indices: 

"3", "9", and "11". By a similar process (repeated application of 
Step 4), the remaining subsets, S^-Sy, can be obtained, and the results 
are summarized as follows: 

S^: 1; S 2 : 2; S^: 3,9,11; S^: 4,10,12; S^: 5; 

Sg: 6,8,14,16; S^: 7,13,15. 

The eigenvalues are evaluated using these submatrices as well as the 
original matrix. The comparative values are given in Table II. 

TABLE II 


Eigenvalues of the K matrix: 


Submatrices 
S^ 0.0 
S^ 0.0 

S- 26.88378000 
41.36168802 
78.46953198 


Original matrix 
0.0 
0.0 

26.88378620 

41.36169439 

78.46953100 
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Fig. 2. 


po 

12 5 6 0 16 7 IJ IS 


Digraph of 16x16 K matrix. 
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11x11 matrix 
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Fig. 3b. Reachability Matrix for Free-Free 
Square Plate after Separation of 
S^,S 2 jS 2 Submatrices. 
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Fig. 3c. Reachability Matrix for Free-Free Square 
Plate after Separation of S^jS^.S^.S^ 
Sub-matrices . 













TABLE II - (Contlnuatxon) - Eigenvalues of the K matrix: 


Submatrxces 

S, 33.24206125 
260.2368175 
333.3811212 
0.0 

S. 12.21269086 
330.7355000 
429.6358695 
1747.485740 

S- 33.24206125 
260.2368175 
333.3811212 

The eigenvalues of the system a 
M=ml, eigenvalues of M“lK=(l/m) 


Original matrix 

33.24206653 

260.2367428 

333.3811150 

0.0 

12.21269085 

330.7355804 

429.6359588 

1747.485839 

33.24206653 

260.2367428 

333.3811150 

e the eigenvalues of M ^ and as 
times eigenvalues of K. 


V. Controllability 


The controllability of a dynamic system represented by X=AX+BU, for 
a given set of actuators can be established using the path matrix ob- 
tained in the eigenvalue analysis. The two conditions to be satisfied 
for controllability are^^: 

(i) All states must be reachable, in a diagraph sense, form 
at least one input. 

(il) The required term rank tests must be satisfied. 


These conditions will be discussed and fully explained in the following 
paragraphs . 


As an example of condition (i) let us consider the system matrices 


A = 


1 1 
0 1 


The digraph corresponding to the A matrix can be represented as 




U 
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Because of the presence of the zero element in the second row of the 
first column, it is observed that the directed graph between the two 
nodal points is not strongly connected. In addition there are also 
directed lines from both nodal points to the points themselves. 

Now consider the control input matrix, B, which has a non-zero 
element in the second row. - i.e. this input directly influences the 
second node and via the digraph has a directed path to the first node. 
Therefore, all states are reachable as a result of the (control) input 
given by B. 

T 

In contrast suppose that the B matrix has the form: B'*' = [1 0]. It 

is then apparent that the second node in the digraph is not reachable 
as a result of this (control) input. 

The term rank [condition (ii)] of a matrix is the maximum rank that 
the matrix may have due to the locations of the non-zero elements. If 
the term rank of the A matrix is equal to the order ^ of the system (n) 
and condition (i) is satisfied , then the system: X=AX+BU is con- 
trollable. If the term rank of the A matrix is less than n, columns of 
the B matrix can be used in place of selected columns of the A matrix in 
an attempt to augment the term rank. The number of columns of the B 
matrix used in this rank augmentation process is an indication of the 
minimum number of actuators required. 

In the example of the free-free square plate the term rank and the 
numerical rank of the K matrix are the same (13) . Accordxng to the 
second condition for controllability, there must be at least three 
actuators to control the plate since the order of the system is sixteen. 

As a particular case the three actuators are selected as shown in 
Fig. A and the controllability is verified. 

The closed loop equations of motion for the system can be repre- 
sented as follows : 



where F = . . . F^g] 

m = mass of the plate 

and the control influence matrix, B, is listed in Table III. System (3) 
is controllable if the pair [-K/m, B/m] is controllable.^^ 
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TABLE III 


Elements of the B Matrix in Eq. (3) 


1.0 

0.866 

-0.198 

-1.1694 

1.2990 

1.125 

-0.2577 

-1.51918 

0.847 

0.73368 

-0.168 

-0.99076 

0.097396 

0.0842 

-0.01932 

-0.1139 


1.0 

0.866 

-0.198 

-1.1694 

-1.2990 

1.125 

0.2577 

-1.51918 

0.847 

0.73368 

-0.168 

-0.99076 

-0.097396 

-0.704 

+0.016129 

+0.095076 


1.0 

- 0.866 

-0,198 

0.117445 

1.2990 

1.125 

-0.2577 

1.52566 

0.847 

-0.73368 

-0.168 

0.99076 

0.097396 

-0.0842 

-0.1932 

0.114387 


For the pair [-K/m, B/m] the reachability condition (i) is very 
easily verified as there is no zero entry in the B/m matrix. The -K/m 
matrix has a term rank equal to numerical rank of 13. By supplementing 
the three columns of the B/m matrix, the -K/m matrix rank is augumented 
to 16 which, together with the reachability condition, assures the con- 
trollability of the free-free square plate for this set of actuators. 


VI. Conclusions 


A graph theory approach for reducing the system matrix to lower 
ordered sub-matrices can be utilized to find the eigenvalues. This 
approach can substantially reduce the numerical effort involved in 
finding the eigenvalues of large order flexible space systems. The 
reachability matrix and term rank concepts are used to verify con- 
trollability and can be computationally more effective than numerical 
rank tests of the system controllability matrix. 
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CHAPTER IV 


THE CONTROL OF LARGE FLEXIBLE ORBITING 
SHALLOW SPHERICAL SHELL STRUCTURES 


The dynamics of shallow spherical shells in orbit are studied in 
Ref. 1. The dumbbell is attached to the spherical shell at the apex to 
achieve gravitational stability of the uncontrolled system.^ In addition 
if it is assumed that the dumbbell is attached to the shell with a spring 
loaded gimballed damping device it can also provide passive restoring 
and dissipative forces. The model developed in Ref. 1 is adopted and 
the effect of the control actuators on the system is incorporated into 
the model. 


INTRODUCTION 

The linearized equations of motion of shallow spherical shells in 
orbit are given by^*: 


4>"+4n^$+(l-n^)iJ;' = 

0"-3n 0-2Ee'lJ’^^£/J = C /J 

ynnl yyyc g 

e"+(fi2-3)e +20 'iJ“Vm Z = SL + -- -- rr 

- I n 1 n M 


(IV-1) 


n c 


n n n 

The linearized equations of motion of the shallow spherical shell 
with the dumbbell in orbit are given by 




N 


<> +4J2 <}i+(l-fj )tf; -c 6 -k 6-E (c e +k e )C^ = C /J 

z'^ z z - znzn y zzc 
n — 1 

Q"-2Q Y+EC^c'4k £ )C^"^= C /J 

y nnlJy y ynTnynz yyc 

I 3X J 

"mV ‘"z”’ 

n n li. 

-(c o'4k 6) -p -2 C^’^^+ E('c e'+k e )C™^+ ECc" e' +v e )C^ = E /M 

z z M^i^y mynynz m znzny nnc 


’'Note: The nomenclature for this Chapter is provided at the end of 

the Chapter. 
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Y"+Cy(l+c^)Y’ + {3+ky(l+c^)}Y+3(l+np9 - (l+c^)E (CyE^+k^e^) C 

+2(Ee'l.^'^^£/l )= 0 
n n 1 y 

6"+7^(1+C2)<S' + {4+k^(l+C2)}5+4(l-fi^)(*> - 

-(l+c^)E(I e’+k e 0 

2 2 n 2 n y 


(n) 


(IV-2) 


The frequencies of the spherical shell are evaluated using the following 
2 

identities. 

y = la [-- g --- (p2-p2)]^ (IV-3) 

2 

where the y's are calculated from 


a r P 
2 ^ 


i „ C ») 






] = 1-v for n= 0,1 


where n represents the number of nodal diameters. 


(IV-4) 


For n>l. Equation (IV-4) must be replaced by a more complex form as 


follows ! 


■4 = 


Sn(u) 

R,(U) 


-1 


(IV-5) 


where 


s„(w) 


4n^ (n^-i) (1-v) {y [J^(y)l^(y)-J^Cp)ij^(y) ] 


+ (n+ 1 ) (1-v) [r (y)- ^ I^(y)] [j;(y)- ^ J^(y)]> 

R^(y) = {(1-v) [yJ^(y)-n2j^(y)]+y2j^(y)H(l-v)n2[yI^(y) 
-I^(y)]-y3l^(y) }-{(l-v)n2[yJ^(y)-J^(y)]+y3j^(y) }{ (1-v)] yl '^(y) 
-n2ln(y)]-y2in(y)} 

Eq. (lV-4) or Eq. (IV-5) is satisfied by an infinite number of the parameter, 
y, for every value of n. For the sample calculations in this chapter. 


we will consider only three such values of y (j = 1,2,) for the cases 
where n = 0,1. 
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For the following numerical values considered (Fig. IV-1) : 

H = Im ; a = 100 m; p = 27.68 Kg/mr^j h = 1cm 

D = Eh3/12(I-v^) 

E = 0.744x10 Kg/mr^; 

R = 5000.0 mrs; = 0.107572254 

p2= 0.252122471xl0"6p4+ 0.107572254 

the frequencies are then calculated, with the results summarized as follows 
Table IV-1 - Frequencies of the Shell (Six Ffodes) 


# of nodal 

# of nodal 



diameters (n) 
0 

circles (j) 
1 

9.076 

0.107593022 

1 

1 

20.52 

0.107678417 

0 

2 

38.52 

0.107946351 

1 

2 

59.86 

0.108475664 

0 

3 

87.82 

0.109516711 

1 

3 

119.0 

0.11114256 


It is observed that all the frequencies are grouped with only a 
difference of 3% between the maximiim and minimum values of the flexible 
frequencies considered in this model. This close grouping is due to 
the dominance of the curvature of the shell on the assumed model. 

The point actuators are modelled as follows. An actuator located 
at (x,y,z) with components (f^jf^.f^) provides the following torques. 


T = yf -zf 

X z y 

T = -xf +zf^ 
y z ^ 

T^ = xf -yf 
z y ^ X 

and the corresponding generic force in the nth mode. 


(IV-6) 
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(IV-7) 


E = / $ • £ dm 
n n 

For the shallow spherical shell it is assiimed that the major elastic 

displacement occurs in a direction normal to the base (y,z) plane. - i.e. 

= <{i^i. Thus for a point actuator located at (x,y,z) 

= ^„(x,y,z)f M (IV-8) 

n n X n 

where <{>^(x,y,z)^ is the n ^ modal shape function evaluated at (x,y,z), 

f represents the component of f in a direction normal to the base plane, 

and M is the modal mass, here considered to be the total mass of the 
n 

shell, since the shell mode shape functions have each been normalized 

with respect to the mass of the shell. 

After defining the state vector for the spherical shell, with and 

without dumbbell, as: 

T 

X = [i}<, ({), 6, e^, e^, e^, e^, e^, ip’ , 6', e^, e^, e^, 

e^, e^, Eg ] (without dumbbell) 

= [y, 6, Tp, <p, 9, z^, e^, y', 6', 'i'*, <!>’, 9', 

e', e’, e! , el, el] (with dumbbell) 
z J 4 0 D 

the equations of motion can be written in the state vector form: 

X = AX + BU (IV-9) 

The B matrix is evaluated using the locations of the actuators together 
with the direciton cosxnes of the actuators . 

Six actuators are assumed to be present and are located on the 
surface of the shell as illustrated in Fig. IV-1. 
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The location of these actuators and the directions of the forces that 


each provides are listed in Table IV-2 



Table 

IV-2 - 


8 

6, 


o 

1 

1) 

90° 

0.28 

2) 

90° 

0.57 

3) 

90° 

0.84 

4) 

270° 

0.28 

5) 

270° 

0.57 

6) 

270° 

0.84 


Actuator Locations and Force Directions 


X 

f. 


y 

0 

0 

0 

0 


0 

0 


^4 

0 


The physical parameters of the shell are the same as those used in the 
frequency calculation. The corresponding non-zero elements of the A and B 
matrices are listed as follows: 


A&B matrices for the case without Dumbbell 


A(I, 1+9) = 1.0 for I = 1,9 

A(ll,2)= 4.0; A(12,3)= 3.0, A(13,4)= -40219.81; A(14,5) = -40251.73; 
A(15,6)= -40351.89; A(16,7)= -40549.78; A(17,8)= -40938.97; 


A(18,9)= -41652.59; A(10,ll)= 1.0; A(ll,10)= -2.0; A(12,13)= 0.02339; 


A(12,14)= -0.0026; A(12,15)= 0.000751; A(13,12)= -0.005846; 
A(14,13)= 0.00066; A(15,14)= -0.00018 


Lower half part of B matrix 


o.onooooE a 
o.noonnn? o 
-n,?oa^atF -•? 

-i 

0.1 -1 

0.a<J<jn;,oF -■« 

OalPTflMF -I 

o.josme -1 

n . t 7< Jfe «6 F -1 


0,onrt0O0F_ o_ 

0.ia7176E -3 

- -0.10717SE 

o.ftoonnnF 0 

0.31?7*1E -3 

0,l660ft'5F 

-O.?a/1T01F -3 

-o.FoaiatE -3 

o.?oa3aiF 

".SnitOTE.-J . 

-0.?3('7tlE„-2 

_ 0.1?055?E 

-n.APJhOfiF -? 

-0,ao6??iE -? 

0.n03«|F 

-0.7AB77nF -? 

0.77a?f,6F -2 

O.poonfcqF 

'’.1 ?'>0S1F _-i 


_-0.1?7f.f>?F 

-o,-^o-?;,qnF -3 

-o.OftOftipe 

•o.?P9l3?f 

-n, 1 /j?o#.QF -t 

o.inn6<J9E »l 

-o.WOfrahE 


-3 

-0.?3a3aoE 


-n,ia7aioE -3 

-3 

0,a6B6«3E 

-*> 

n.3t325?E’-3 

-3 

0.7353'5f.E 

• a 

-O.tPaaiiE .a 

.-L . 

0.5*16<J7E 

-? 

-0,?3«7ilE -? 

-1 

-0,fr?aftooE 

-? 

-n,a9^,3?iF -? 

-3 

-0.7^«770E 

-2 

n,77a?6(>E -? 


,-0.t3Pa':nE 

-1 

_^,loqja?F -? 

-1 

O,3036O0E 

-3 

o.^hOftinE -? 

-1 

0. toaofeOF 

-1 

-0.ino<.99E -1 
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Non-zero elements of A with dumbbell 
A(I, I+ll)= 1.0 for I = 1, 11 

A(12,l)= -193.0; A(12,5)= 6.0; A(12,7)= 576.32; A(12,9)= 3647.6; 

A(12,ll)= 6088.63, A(12,12)= -2.7568, A(12,17)= -0.0234; 

A(12,18)= 8.365; A(12,19)= -0.000751; A(12,20)= 52.93; A(12,22)= 88.3432; 
A(13,2)= -194.0; A(13,3)= 2.0; A(13,4)= -8.0, A(13,13)= - 2. 757; 

A(14,15)= 1.0; A(15,2)=100.0; A(15,4)= 4.0; A(15,13)=1.451; 

A(15,14)= -2.0, A(16,l)= 100.0; A(16,5)= 3.0; A(16,7)= -303.33; 

A(16,9)= -1919.8; A(16,ll)= -3204.54; A(16,16)= 1.451; A(16,17)= 0.0234; 
A(16,18)= -4.4038; A(16,19)= 0.000751; A(16,20)= -27.8554; A(16,22)= -46.496; 
A(17,6)= -40216.81; A(17,16)= -0.00585; A(18,l)= 75.834; A(18,7)= -56712.312; 
A(18,12)= 1.10032; A(18,16)= 0.00066; A(18,18)= -238.88; A(19,8)= -40348.9; 
A(19,16)= -0.0001878; A(20,l)= 479.97; A(20,9)= -144747.31; A(20,12)= 6.9641; 
A(20,20)= -1511.9; A(21,10)= -40935.97; A(22,l)= 801.16; A(22,ll)= -215581.68; 
A(22,12)= 11.6245; A(22,22)= -2523.672. 

The B matrix of the shell with the dumbbell is 

the lower part of the B matrix of the shell without the dumbbell. The locations 

of the non-zero elements in the matrices of the spherical shell with and 

without the dumbbell, respectively, are shown in Figs. IV-2 and IV-3. 

4 

Control laws are selected using pole clustering techniques for the 
spherical shell without the dumbbell. This same control law is then applied 
to the case of the shell with the dumbbell, i.e.- it is assumed that the 
dumbbell position and rate information .may not- be directly observable and 
is not included in the control law. (As long as' the dumbbell-shell system 
is controllable, and stabilizable, the dumbbell will return to its desired 
local vertical equilibrium orientation after the transients have been 
removed . 
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It is possible to design a control law for the dumbbell-shell system 
that results in a controllable-stabilizable system without including 
the dumbbell information within the control law.) 

A typical time history of the required control forces for the shell 
without the dumbbell is shown in Fig. IV-4 for non-dimensional initial 
position displacements in all state components of 0.01, with the control 
law based on placement of all the poles so as to have a dimensionless 
negative real part of (-1.72). This will provide a system time constant 
of 461 seconds in all modes. A similar response was generated using 
the same control law for the case of the shell with the dumbbell where 
the dumbbell inertia ratios (c^=C 2 ) were assumed to be 0.9, the dumbbell 
spring constants, k =k =100, and damping coefficients selected so as to 

y z 

provide 0.1 critical damping. A comparison of the maximum force amplitude 
in each of the actuators and also of the total force impulse required 
is given in Table IV-3 for the two cases. It can be seen that although 
there is little difference in the maximum force amplitudes, there is 
approximately a 25 percent savings in fuel consumption by using the shell- 
dumbbell system. 
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Table IV 3 ; Maximum Force Amplitudes of the Shell with 6 Actuators 


-p- 

00 



Without 

With 



dumbbell 

dumbbell 


^1 

103.411 

103.410 

Initial Conditions 

^2 

105.398 

105.398 

*0 ■ *0 ” *'0 “ ‘l(0) " ■ 1=6(0) ■ 

^3 

317.83 

317.83 

«0 ■ *0 ■ *0 ■ 'l(0)= • • • ==6(0) ■ O-O 

^4 

107.55 

107.55 

with dumbbell 


77.09 

77.09 

^0 = ^0 = ^0 = ^0 = 

^6 

127.32 

127.318 

( 


Total 

Force impulse: 432.8 N-Seconds 320.49 N-Seconds 25% saving with dumbbell. 
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H = lin h = 1cm a = 100 m 
m = 869679.8058 kgs. 




NOMENCLATURE - CHAPTER IV 


: modal amplitudes 

: external torques with components 

: co-efficients of viscous damping 

: c/J (jJ , c /J (0 ,c /J 0) , respectively 

y c’ y y c z z cJ ^ 


( 0 , 0 ), 


( 0 , 0 ), respectively 


c(m)c(n) 

y y 


c(m)c(n) 

z z 


J /I, 

y d 

Z d 

J /m£^ 

y 

modal component of external forces 


moment of inertia of the dumbbell 


principal moments of inertia of the undeformed shell 

torsional spring constants 

k/J u)^, k /J oj^, k /J 0)^, respectively 
ycyyczzc 

characteristic length 
length of the dumbbell 
modal mass of n^^ mode 


mass of shell 


tip mass of the dumbbell 


radius of curvature of the shell 


: vector from the mass center of the shell to the 

origin of (x^,y^2^) 

: position vector of a generic point on undeformed 

shell (x y z ) 
c c c 

: radial distance of the shell element from the 

symmetry axis of the shell 
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t 


e 

(Y,5) 


e 

n 


e,4>,<j) 

p 

T 

T 



0) 


01 

c 


0) 

n 


Q , Q 
X y z 


(•) 


( )’ 


time 

polar angle 
phase angle 

dumbbell deflection angles 


A /Z 
n 

damping ratio 

pitch, yaw and roll angles, respectively 

mass density 

CO t 
c 

coordinate frame, (x y z ) 

c c c 

dumbbell frame (x,yjZ.) 

dad 


body angular velocity vector, 
orbit angular velocity 
natural frequency of n*"^ mode 


(o) 0) 0) ) or (oi (0-0) ) 
X y z r e X 


: (J^-Jy)/J^, ’ respectively 

: d 

dt 

d 

' dT 
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CHAPTER V 


THE DYNAMICS OF LARGE FLEXIBLE 
EARTH POINTING STRUCTURES WITH A 
HYBRID CONTROL SYSTEM 


INTRODUCTION 

Proposed future applications of large space structures include; space 
based power generation and transmission (to earth) ; communications ; 
earth resource observation missions; and electronic mail systems. The 
applications described here all reqxiire that the largest surface (length) 
of the system be nominally oriented along the orbital tangent or normal 
to the local vertical. In order to satisfy the mission requirements for 
some of these applications control of the shape of the structure as well 
as the orientation will be required. 

Previously, the equations of motion of a general arbitrary flexible 
spacecraft in orbit have been developed. ^ As a specific example, the 
dynamics and stability of a long, flexible beam constrained to move 
only in the orbital plane without control was considered, with the prin- 
cipal emphasis placed on the motion about the nominal earth pointing 
(local vertical) orientation.^ 
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It was observed that for small amplitude pitch and flexural oscillations, 
the rxgid pitch mode was not affected by the elastic modes, and that 
the elastic motion was coupled to the pitch mode and described by a 
Mathieu equation. The possibility of parametric instability at very 
low natural elastic frequencies was demonstrated. 

In a recent paper^ the work of Ref. 2 was extended to consider: 

(1) the motion and stability of the beam about a nominal local hori- 
zontal orientation; and (2) the dynamics of free-free homogeneous plat- 
forms in orbit with emphasis placed on a circular plate structure. Both 
models included only the presence of an articulated passive damping de- 
vice. 

It is well known that, due to an unfavorable moment of inertia 
ratio, a rigid orbiting thin beam, when displaced from a local hori- 
zontal orientation, will not tend to return to this orientation in the 
absence of external control forces. This tendency is also observed 
when the flexural dynamics are added to the model. Such a structure 
can be gravitationally stabilized by using a rigid dumbbell with 
proper moment of inertia (Fig. 1). In Ref. 3 the dumbbell was assumed 
to be attached at the center of beam by a spring loaded hinge with 
viscous damping also present. It was seen that the pitch motion, the 
dumbbell librational motion, and the anti-symmetric elastic modes of 
the beam were all coupled. Within the linear range, the beam's symme- 
tric elastic modes were decoupled from the pitch and dumbbell motions. 
The frequencies of the symmetric elastic modes were adjusted to include 
the effect of the dumbbell. Using the Routh-Hurwitz criteria, stability 
* conditions were derived and limits on key system parameters established. 
A parametric study of the least damped mode characteristics for dif- 
ferent system parameters was carried out. In Ref. 3 the three dimen- 
sional equations of motion for a flat plat" in orbit with a two-degree- 
of- freedom gimballed dumbbell were also developed and, within the li- 
near range, it was seen that only certain of the elastic modes were 
influenced by the dumbbell motion. Thus, in order to damp the motion of 
such systems nominally oriented along the local horizontal in all the 
modes generally contained within any trxmcated system model, it was 
concluded that the use of active devices would also be required. How- 
ever, it was thought that with the use of the passive gimballed dumb- 
bell together with the active controllers, the fuel consumption could 
be significctuuly l educed without exceecing tne previous limits on trie 
peak forces. 

The present paper examines the dynamics and control of such large 
flexible earth pointing structures in orbit \mder the influence of a 
hybrid control system consisting of the passive gimballed dumbbell to- 
gether with active point thrusters. 
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A THIN UNIFORM BEAM IN ORBIT NOMINALLY ORIENTED ALONG THE LOCAL 
HORIZONTAL 

The Mathematical Modelling of a Flexible Orbiting Beam - Local 
Horizontal Orientation 


The equations of motion for a thin uniform beam in orbit with its axis 
nominally along the local horizontal (Fig. la) is developed in Ref. 3. 
The beam is assumed to undergo only inplane angular motions and defor- 
mations and it is assumed also that the center of mass of the beam 
follows a circular orbit. The beam's elastic motions are considered 
to be unconstrained and the longitudinal vibrations of the beam are 
assumed to be negligible in comparison with the transverse vibrations. 
For the case of small amplitude pitch oscillations of the beam, the 
linearized equations of motion are derived as^. 


where 


0"-36 = C /J (J^ 
y y c 

e" + [({ij /o) )^-3]e = E /M 


n 


n 


n c 


( 1 ) 

( 2 ) 



mdeformed beam length 
pitch angle 

orbital angular velocity 
external torque about the pitch axis 
beam pitch axis moment of inertia 
n*"^ flexural modal amplitude 

A^/2. nondimensionalized n^^ flexural modal amplitude 

n^^ modal mass 

where x = m t 
dx c 

time 


As can be seen from Eq. (1), the uncontrolled local horizontal orien- 
tation of the beam represents an imstable motion. This unstable con- 
figuration of the beam can be stabilized either by providing active 
control torques or by adjusting the moment of inertia distributxon of 
the system such that the resulting gravity-gradient torques provide 
stabilization. A gravitationally stabilized system using a connected 
rigid dumbbell will be considered in the next section. Eqs. (1) and 
(2) can be represented in a state vector form for a ready application 
of linear control theory. 
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After defining the coordinates according to: 


e = x^, 

"l = ^2’ 

^2 ^3 * • • • • j 


(3) 

and e ' = X , . , 
n+1’ 

" \+2’ 

£ * “X 

2 n+3^ * 

“ =‘2n’ 



Eqs. (1) and (2) take the form 


X = AX + B U 
c c 


where X 


* [x^, X, 


j ♦ • • 


n+1 


"'Zn^ 


(4) 


A = 


where 


0 

'V; 


0 

'\i 

B = 

c 

U = 


0 


nxn null matrix; = nxn identity matrix 

[”3 -(uj-/u) )^+3 .... -(oj -/(ij )^+3] 

1 c n-1 c ■' 

0 = null vector of order n 

1 = identity vector of order n 

[Cy/Y^, Ej_/M^032i, . 


0 

.IJ 


(5) 


..2 T 7 /M ...2 0 _ £ 

n-1 n-1 C-* 


Dtrabbell Stabilized Flexible Beam in Orbit with the Longitudinal 
Beam Axis Nominally Along the Local Horizontal 

A passive stabilization of the beam can be obtained by using a 
rigid dumbbell with proper moment of inertia. In Ref. 3, the equa- 
tions of motion for a beam with a dumbbell assumed to be attached at 
the center of mass of the beam (Fig. lb) through a spring loaded 
hinge and halting vj.scous rotational damping have been developed. In 
addition to the assumptions made in developing Eqs. (1) and (2), it 
is further assumed that the dinnbbell mass is concentrated at the 
tips and that the viscous force at the hinge is linear. With the 
usual assumptions of small pitch amplitude and dumbbell oscillations 
and flexural deformations, the linearized equations of motion in 
the absence of active control and external forces are obtained as^: 

e"-t^9'+(k-3)6-'ca’-ka+Z(«'+ke = 0 

n n n z 

a"+c 7a'+(c k+3)a-c^'ce'-c^ke-I (cc’+ke = 0 (7) 

^ J- innn lz 
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= 0 (8) 


e"+(n^-3)e - {k(a-0)+c(a’-e’)}C^’^^J /M £^+I(ce'+ke 
nnn zynminmz 

where = J (m,n = 1,2,...) and M = mass of the 

z ^ 2 z n XI 

beam for all n. 


k = k/J u)^; c = c/J 0 ) 
y c y c 

k = torsional restoring spring constant at the hinge 


c 

a 


X 


= viscoxis damping coefficient 


.(n) _ 


= angle between the dumbbell axis and the local vertical 


3z 


2=0 

= beam shape function 


c^ = J /l^, = pitch moment of inertia of the dumbbell. 


= u /lO 

n c 


The following observations can be made from a study of Eqs. (6) -(8) : 

(a) the pitch motion of the beam, the dumbbell motion (o) , and the 
elastic motion of the beam (e ) are all coupled to each, other; (b) 
within the linear range the elastic modes for which - 0 (the 
symmetric modes) , are completely independent of the pitch and dxlmbbell 
motions, furthermore, these modes do not influence either the pitch 
or dinnbbell motion; (c) becuase of the presence of the dumbbell the 
natural frequencies and the mode shapes of the symmetric modes of the 
beam differ from those of the free-free beam. The frequencies and mode 
shapes can be easily obtained by replacing the dumbbell by a concen- 
trated mass, equal to that of the dtimbbell, at the center of the beam.** 
Eqs. (6), (7) and (8) can be rewritten in the state vector form of Eq. 
(4). With the inclusion of the first two free-free beam elastic modes, 
the eighth order matrix A has the form 


A 


where 

0 



^ ^2j 

4x4 null matrix 


4x4 identity matrix 





Modelling of the Point Actuators for the Beam 


Let there be p actuators located along the beam at distances 
K from the center of mass of the beam and the actuator forces be f^, 
f^,...f , respectively. For small elastic displacements, the control 
torque Hue to the actuator may be expressed as 



( 10 ) 


where 

f = k f 6 (5-5.) ; i,j,k 

J J J 

r - 5ji + qj(lj) 

q. = ? ♦" (5) A^(t) k 

^ n=l ^ “ 


itnit vectors along the 5>n»5 
directions. 


n 

and (j) in the z component of the modal shape function corresponding 
to the n^“ mode. 
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For small deflections 0, and there are no force components along 
the X and y directions.^ Therefore, 

\ - 

and for p actuators 

iip = I/Pj = -J[flCi+f2l2+....+fpSp] - (12) 

The control forces can be transformed into the corresponding modal 
forces by^»^ 

Enj = / ? "" • fjdm (13) 

under the assumptions previously stated 

E^j •= fj <}i^ (5j) X (a constant) (14) 

The constant can be incorporated into f, so that the effect of all p 
thrusters on the n^^ generic mode can b^ expressed by 

( 15 ) 

The control vector, U , can now be related to the actuator forces, 
actxiator locations , and modal shape functions by 


U = 
c 


<t.l(Ci)/KL£aj2, <{,1(C2)/M3^^2^ 



I 





» 


nxp 

Control of the Orbiting Beam System 


pc 

)/M .iUjj2 

n-1 c 



pxl 

( 16 ) 


Since the beam with its axis nominally alone the local horizontal 
represents an unstable motion, active control using thrusters will 
have to be provided. In the case of the beam stabilized using the 
dumbbell, active control is still required to achieve an acceptable 
time response during attitude corrections and to provide damping 
directly into the symmetric modes. 


A single actuator located at the left end of the beam is assumed to 
provide this control. The fimdamental natural frequency of the free- 
free beam is taken as 0.01 Hz. and the center of mass is assumed to 
move in a 250 nautical mile altitude circular orbit. For this case, 



As an example, a 100 m. long slender hollow tubular beam of aluminum 
(2014T6) with outside diameter of 10.79 cm, and thickness 1.06 cm 
would exhibit these frequencies. 

The controllability of the beam without the dumbbell Eq. (1) has been 
studied in Ref. 5. It is shown that for a beam with an actuator at 
one end the necessary and sufficent conditions for controllability 
are (with two modes in the model) 

0)^ 5^ (1^2 ; /~3 (t>^ ; ^ 

The controllability of the beam with the dumbbell was verified using 
the ORACLS® computer algorithm. 

First, the control law is selected for the case of the beam without 
the dtimbbell. Then the gain values are determined so that the system 
response is acceptable. The same control law and gain values are then 
applied to the case of the beam with the dumbbell and the system re- 
sponse is obtained for comparison. 

The rate feed-back of all the states is required in order to obtain 
damping in the system. Further, as can be seen from Eq. (1), there 
is a pole on the positive real axis indicating instability. Hence, to 
move this pole to the left of the imaginary axis of the complex plane 
a pitch position feed-back is necessary. Assuming that pitch, pitch 
rate and modal amplitude rates are fed-back, the closed-loop system 
characteristic equation for the beam with the two modes included in 
the model is derived as 

s®- {aK^+b(K5+Kg)s5 + {c+d-CaK^^+S) }s‘* 

- {bCdK^+cKg) + aK^(c+d)-3b(K2-HCg)>s3 

+ {cd-(aKj^+3) (c+d) }s^ - {acdK^-3b(dK2+cKg) }s 

- (aK^+3)cd = 0 (17) 

where 

a = i'/ 2 JyU )2 ; b = 

c = (o) /lo )2- 3 ; d = (io„/o) )^- 3 
1 c z c 

and K, = pitch gain; K,,K,-,K- rate gains for pitch and modal 
amplitudes 

By inspecting the coefficients of Eq. (17) , the necessary conditions 
for stability of the system are found to reduce to the following, 
noting that a,b,c, and d are all positve for the beam vmder consider- 
ation, 

(i) aK^+3 < 0 or < -3/a 

(ii) K^, K^, Kg should be negative. 
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Numerical Examples - Orbiting Beam System 


By selecting = -0.07, = -0.005, = 0.1 and K, = -0.03, the 

system response was obtained using numerxcal integration as shown 
in Fig. 2. Initial errors of 5° in pitch and 1.5m in the two modal 
amplitudes were assumed. The pitch motion and the two modal oscil- 
lations are seen to be under control and the transient motion has 
been essentially removed in about 6 orbital periods. The maximum 
control force required is about 0.025 Newtons. 

The same control law with identical gain values is used to obtain 
the response of the dumbbell stabilized system. No feed-back from 
the dumbbell is assumed for the sake of simplicity. The following 
dumbbell parameters are assumed: 

Jy/I^ =0.9 ; k = 60 ; 7 = 0.562. 

Fig. 3 shows the system response obtained by numerical integration for 
the dumbbell stabilized system with the same initial conditions as in 
Fig. 2. In comparsion with Fig. 2, the pitch and the second modal 
oscillations of the dumbbell stabilized beam are seen to have better 
response characteristics. The dumbbell oscillations directly affect 
the asymmetric mode (second mode) and as a result the as 3 nnmetric 
modal oscillations are very effectively damped. The symmetric mode 
oscillation (first mode) is not affected by the presence of the dumb- 
bell and is controlled mainly by the actxiator. The pitch motion is 
also influenced by the dumbbell motion as long as there is relative 
motion between the dumbbell and the pitch mode - i.e. (6-a) is not 
zero. The peak control force required is about 0.012 Newtions. There 
is a total savings of about 45 percent in the force-impulse require- 
ment as compared to the case of the beam without the dumbbell. The 
change in the frequencies of the beam due to the addition of the 
dumbbell was determined using the results of Ref» 4. Using these 
frequency values the response of the system for the same initial con- 
ditions was obtained. No significant change was seen from the results 
in Fig. 3. Other numerical results were obtained with different ini- 
tial conditions (not shown here) indicating , in each case, better re- 
sponse and fuel savings for the dumbbell stabilized system. 

A THIN UNIFORM PLATE IN ORBIT WITH ITS MAJOR SURFACE NOIUNALLT IN IHE 
LOCAL HORIZONTAL PLANE 

A flat plate with its major surface in the local horizontal plane is 
gravitatiionally unstable in the absence of external restoring torques. 
In order to provide stabilizing torques, active control using thrusters 
can be employed or the moment of inertia distribution can be readjusted 
using a rigid dtmbbell, as in the case of the beam. A rigid dumbbell 
is assumed to be attached at the center of mass of the plate (Fig. 4) 
by a spring loaded double gimballed joint. Here, the dumbbell would 
have two degrees of freedom as described by the angles y and 6 as 
shown in Fig. 4. 
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The equations of motion for a thin uniform flat plate with its major 
surface nominally in the local horizontal plane and carrying a rxgid 
dumbbell at the center of mass have been derived in Ref. 3. The assump- 
tions made in the derivation can be summarized as: (a) the plate is 

thin and has uniform mass and stiffness properties; (b) the center of 
mass of the plate is moving in a circular orbit; (c) the longitudinal 
vibrations of the plate are negligible; (d) the plate is uncon- 
strained; (e) the mass of the dumbbell is concentrated at the tips of 
the dumbbell; (f) the center of mass of the plate and that of the 
dumbbell coincide; (g) the attitude angles and vibration amplitudes 
are small. With the above assumptions one can derive the linearized 
equations of motion in the following nondimens ional form^ : 


r' - n ii»-(l+Q^)<j)* = 0 

X & 

<J>" - 6'-k 6+2 (c e’+k e = 0 

2 zn znzn y 

e'« - 30-7 Y+ZU t'+k e = 0 

y y n y n y n z 

e" + (n2-3)e -{(c y'-Hc y)J-(j‘’^V(7 6’+k 6) J 


y z 


z y 


+2{(c e*+k E )C^““^+(c e'+k e = 0 

m ymymz zmzmy 


(18) 

(19) 

( 20 ) 

( 21 ) 


y" + c„(l+c,)Y'+{3+k (14c,) }Y+6e-(l+cj2(c, E'+k E^)cf“^= 0 (22) 


I'n' y n y n' z 

5" + c^(l4c2)6’+{A+k^(l+C2)}6+8<j)-2i|;’-(l+C2)Z(I^e^ 

+k E 0 

z n y 


(23) 


where, \p, <J>,. 6 = yaw, roll and pitch angles, respectively; e = 

= modal amplitude; M = modal mass; £ = characteristic length 
(eg. radius of a circular plate or length of the side of a square 
plate); y, 6 = dximbbell angles; J ,J ,J = principal moments of 
inertia of the plate; = principal^moment of inertia of the dumbbell; 

= (J -J )/J ; k ,k = torsional spring stiffness; c ,c = damping 
co-efficients ;^and^ ^ ^ 
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03^ = n^^ natural frequency of the plate 

03^ = orbital angular velocity 
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It is interesting to note that the elastic modes for which = 0 

and ~ decoupled from the ip, <f>, 6, y and 6 motions. 

The Case of a Flexible Square Plate 

As a specific example of a flat plate, a thin homogeneous square 
plate is considered. The natural frequencies and mode shapes of a 
free-free square plate have been determined in Ref. 7. The parameters 
of the plate were assumed to be as follows : 

material of the plate is aluminum; sides of the plate = 100m; thick- 
ness of the plate = 0.01m; density = 2768.0 kg/m^; Young's modulus = 
0.7441x10^° kg/m^; Poisson's ratio =0.33. 

The first five nodal patterns along with the frequencies obtained 
for the above plate parameters are shown in Fig. 5. It may be noted 
that the slopes at the center of the plate in all mutually orthogonal 
directions for the first three modes are zero and therefore these 
modes are completely decoupled from the rigid body rotations 
and the dumbbell motions within the linear range. Hence to see the 
effect of the dumbbell on the dyn atm* cs of the plate it is essential 
to include at least the fourth and fifth modes of the plate. 


Control of the Flexible Orbiting Plate 

In order to control the doubly symmetric modes (the first three modes 
in the present model) and to obtain better time responses in 
and the as 3 nnmetric modes, active control using thrusters will be em- 
ployed. Since the plate is homogenous with sides of equal length, 

J = J and hence =0, resulting in an open-loop pole at the origin. 
T^is problem will be addressed, in this paper, by providing a separate 
actuator, f^_ (Fig. 4) which can control the yaw motion independently, 
in addition to an actuator which can control the other degrees of 
freedom, £2 (Figs. 4 and 5). It is assumed that actuator, f^, has 
a thrust direction nominally in the plane of the undeflected plate, 
whereas actuator, f 2 » has a thrust direction normal to the plane of 
the undeflected plate. f 2 is selected so that it does not lie on any 
of the nodal lines as shown in Fig. 5. 


The effect of the actuators can be modelled as follows; 


Let t = f_l + t 1 + 

! y ( 24 ) 

r = 5i + rij + ?k 

AAA 

where r defines the postion of the actuator on the plate and i, j, k 
are unit vectors along the principal axes of the plate. Then, 


T = rxf (25) 

The generic force due to the p actuator on the r mode is given by® 


= // W^(S,n)k • 6(5-5p,n-Tip)fp(t)kd5dn 


(26) 


th 


where W^(5p,Tip) is the r modal shape function of the deformed plate. 
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It can be seen that both actuators, and £ 2 , will produce torques that 
affect the rigid modes, Eq. (25), but that only £2 could produce a ge- 
neric force on the flexible modes, Eq. (26). For this applicatxon 

Er = W^(523h2)f2 (27) 

where ho are the coordinates of f. in the plane of the undeflected 
plate. 

The equations of motion can be rewritten in the state vector form as in 
Eq. (4). The components of the state vector can be identified as, 

= 41, X2 = <}>i = 6 , x^ = Y. x^ = 6 , Xg = e^, = ^ 2 , 

^8 ^ ^ 4 ’ ^9 ^ ^ 4 ’ ^10 “ ^ 5 ’ ^11 ^ ^20 ^ 5 " 

The corresponding matrix, A, of order 20x20 is obtained in the form of 
Eq. (9). For the selected location and number of actuators the matrix 
B^jOf Eq. (4) becomes 




, 0.96 0 0 

^1“ 

L 0 1.36 -1. 


09 0 0 0.019 0.011 0.018 0.013 0.0131 


Parameters of the dumbbell are taken as c^^ = C 2 = 0.9; k • = k =60.0 
and c” = "c =1.12. For this system the controllabilit 3 ^was established 
by using tfie ORACLS® computer alogorithm and the control U = KX is se- 
lected such that (A-BK) has the identical negative real part in each 
of its eigenvalues. This algorithm was found to be useful in deter- 
mining the stabilizing gains of the system. The gains were first deter- 
mined for the case of the plate without the dumbbell. These gain values 
were then used to control the dumbbell stabilized plate. No feed-back 
from the dumbbell motion was included. For initial conditions of 5° 
in pitch and 1.5m in each of the modal amplitudes, time responses were 




the dumbbell. 


Fig. 6 shows the time response of the square plate without the dumbbell 
but under active control. After an initial overshoot in pitch, roll 
and yaw, respectively, the attitude errors are seen to approach zero 
is about six orbital periods. The first three modal amplitudes are 
damped out in about two orbits. The fourth and fifth modsd amplitudes 
require about 6 orbital periods to be completely removed. The 
maximum control force required for the yaw actuator, fj_, is about 3 
Newtons and that for the pitch-roll actuator, £ 2 , is about 30 Newtons. 
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In Fig. 6 only the combination of the absolute values of fj and £2 
depicted. The total combined force impulse required is approximately 
48000 N-sec. 

Fig. 7 shows the time response for the dumbbell stabilized plate with 
active control and with the same gain values and initial conditions for 
the case of Fig. 6. The time response is shown for only two orbital 
periods in which the pitch, roll and yaw motions have been damped out 
effectively. The significant improvement in response of the fourth 
and fifth modal amplitudes as compared with Fig. 6 should be noted. 

These two modal responses (asymmetric modes) are controlled mainly 
by the dumbbell motion. The maximum control force for the yaw actuator, 
f]_, is about 3 Newtons and that for the pitch-roll actuator, f 2 > is 
about 25 Newtons. In addition to the improvement in the responses, 
the dumbbell stabilized plate requires about 15 percent less total force 
impulse than that for the case of the plate without the dumbbell (Fig. 6) 
for the initial conditions and the control law considered here. 

CONCLUSIONS 

The dynamics and control of large flexible structures nominally oriented 
along the local horizontal are analysed. A thin, long flexible beam 
and a thin flexible square plate are considered for specific appli- 
cations. Orientation and shape control is assmed to be provided by 
a hybrid system consisting of a passive spring loaded gimballed dumb- 
bell damper in addition to active point thrusters. The performance of 
the hybrid controlled structures is compared with that of the struc- 
tures having only active control. In general, substantial Improvement 
in response and fuel consijmption for the case of the hybrid systems 
is shown. 
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Fig. la. Uniform Flexible Beam Nominally Oriented Along the Local 
Horizontal. 


X 



Fig. lb. Dumbbell Stabilxzed Flexible Beam Nominally Oriented 

Along Local Horizontal wxch Passive and Active Controllers. 
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Fig. 2. Time Response of Flexible Beam with Active Control and 
No Dumbbell. 
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Fig. 3. Txme Response of Dumbbell Stabilized Flexible Beam 
with Active Control, 
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Fig. 6. Time Response of Flexible Sqviare Plate with 
Active Control and No Dumbbell. 
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Fig. 7. Time Response of Dumbbell Stabilized 

Flexible Square Plate with Active Control. 
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VI. GENERAL CONCLUSIONS AND RECOMMENDATIONS 


The effect of solar radiation pressure as the principal environ- 
mental disturbance torque has been incorporated into the model of the 
orbiting rigid shallow shell. For low earth orbit applications, computer 
simulation results indicate that within the linear range the rigid modal 
amplitudes are excited in direct proportion to the area to mass ratio. 

These amplitudes can be further amplified if the system is designed 
such that any of the natural frequencies are at the fundamental (or basic 
harmonics) of the forcing frequencies - i.e., orbital angular velocity 
or twice orbital angular velocity. 

Graph theoretic techniques can be used for reducing the system 
(stiffness) matrix to lower ordered submatrices which can facilitate 
the calculation of eigenvalues for large ordered systems typical of 
LSST structures . The related reachability matrix and term rank concepts 
can be used to examine controllability with a reduction in the computations 
normally required to ascertain the rank of the system controllability 
matrix. 

The dynamics and control of large flexible beams, plates and shells 
nominally oriented along the local horizontal are analyzed. Orientation 
and shape control is assumed to result from a hybrid system comprised of 
a passive spring loaded gimballed damper in addition to active point 

thrusters. In general, substantial improvement in response characteristics 
and fuel consumption is realized for the properly designed hybrid system 
when compared with the active system operating alone. 
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It should also be emphasized that the passive dumbbell system is reusable 
so that the savings in fuel consumption represented by a single transient 
response would be multiplied many times during the total mission lifetime. 

Possible extensions to the work reported here could be: (1) further 

analysis of environmental effects on LSST systems to include the inter- 
action of the environmental torques with the effects of flexibility; (2) a 
study of the sensor and actuator dynamics to include actual expected delays 
in these systems together with noise characteristics expected to be present 
in the plant as well as the measurement devices; (3) further application 
of the graph theory approach to consider systems with damping present and 
to aid in the calculation of eigenvectors as well as the eigenvalues pre- 
viously considered; and (4) the adaptation of the previously developed 
general model of a large orbiting flexible system to more complicated, 
non-lstropic systems such as the proposed Hoop-column Maypole configuration. 


6.2 



End of Document 



